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Abstract. We obtain strong and uniform asymptotics in every domain of the complex 
plane for the scaled polynomials a(3nz), b(3nz), and c(3nz) where a, b, and c are the type 
II Hcrmite-Pade approximants to the exponential function of respective degrees 2n + 2, 
2n and 2n, defined by a{z)e' z - b(z) = 0(z 3n+2 ) and a(z)e z - c(z) = 0(z 3n+2 ) as 
z — > 0. Our analysis relies on a characterization of these polynomials in terms of a 3 x 3 
matrix Riemann-Hilbert problem which, as a consequence of the famous Mahler relations, 
corresponds by a simple transformation to a similar Riemann-Hilbert problem for type 
I Hcrmite-Pade approximants. Due to this relation, the study that was performed in 
previous work, based on the Deift-Zhou steepest descent method for Riemann-Hilbert 
problems, can be reused to establish our present results. 



1. Hermite-Pade approximation 

In this paper we consider quadratic Hermite-Pade approximation to the exponential 
function. Type I quadratic Hermite-Pade approximation to the exponential function near 
consists of finding polynomials Pm,n 2 ,n 3 , Qni,m,n 3 and r ni n2 n3 of degrees rii,n 2 and n 3 
respectively, such that 

Pn-\,n<2,n;\ (^) ~\~ Qm,n2,n 3 (^) ^"n\,n2,n 3 (^0^ G{z ' ), Z ► 0. 

If we set the right hand side equal to zero and solve for e z , then we obtain the algebraic 
function 

— Q , ni,n2,«3(^) ^ \j Qni,n2,n 3 ~ ^Pni,n2,n 3 { z ) r n 1 ,n2,n: i { z ) 
2?"ni,n2,n3 (^) 

as an approximation to e z . Type II Hermite-Pade approximation is simultaneous rational 
approximation to e~ z and e z and consists of finding polynomials a ni>n2)n3 , b ni>n2;ri3 and 
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c ni ,n 2 ,n 3 of degrees at most n 2 + n 3 + 2, ri\ + n 3 and ri\ + n 2 respectively, such that 

^ni ,n,2,rt3 (^)^ ^ni-JiiSki (^0 0(^Z ), Z > 0, 

Qjll,na,li3 (^0*^ ^ni,ri2,W3 (^) (D{z ), £ > 0. 

This gives the rational approximants ^,713^3 (^/ani.na.nsC^) to e" 2 and ^^^(^/^^^^(z) 
to e z , and both rational approximants have the same denominator. It is well-known that 
for the case of exponentials, all indices ni, n 2 , n 3 are normal, i.e., the polynomials a ni ,n 2 ,n 3 , 
T12 n>3 and c ni ,n2,^3 exist and are unique up to a normalization constant with degrees 
exactly n 2 + n 3 + 2, ni + n 3 and n\ + n 2 , see jTTJ Theorem 2.1, p. 129]. 

Hermite-Pade approximation to the exponential function have been of interest since 
Hermite and have recently been investigated by Borwein jl] , Chudnovsky |5 a , and Wielon- 
sky [22 ■ The asymptotic distribution of the zeros for the scaled type I Hermite-Pade 
polynomials, 

P n {z) = Pn,n,n(3nz), Q n {z) = q n>n) n{3nz), R n {z) = r n ^ n (3nz) (1.2) 

and their asymptotic behavior as n — > 00, have recently been studied in detail by Stahl (T^l 
[TTHIHIHH] and by Kuijlaars, Van Assche, and Wielonsky j^j, see also jHJ. In a Riemann- 
Hilbert problem for type I Hermite-Pade approximation was formulated. The asymptotic 
analysis of this Riemann-Hilbert problem with the Deift-Zhou J7| steepest descent method 
for oscillatory Riemann-Hilbert problems and Stahl's geometric description of the problem, 
allowed the authors of [H] to find strong asymptotic formulas for the polynomials (jl.2j) as 
well as for the type I remainder term 

E n (z) = P n (z)e- 3nz + Q n (z) + R n (z)e 3nz (1.3) 

that hold uniformly in every region of the complex plane. The paper [S] contained the first 
instance of a steepest descent analysis for a 3 x 3 matrix valued Riemann-Hilbert problem. 
It was followed by the papers PjJ HI] which dealt the asymptotic analysis of 3 x 3 matrix 
valued Riemann-Hilbert problems arising in random matrix theory. A Riemann-Hilbert 
analysis for rational interpolants for the exponential function was carried out in |24j . 

It is the aim of this paper to show that the Riemann-Hilbert analysis of [5] also produces 
the corresponding asymptotic results for the scaled type II Hermite-Pade polynomials 

A n (z) = a n ^ n (3nz), B n (z) = b n>n>n (3nz), C n (z) = Cn tn>n (3nz), (1.4) 

and for the type II remainder terms 

EW(z) = A n (z)e- 3nz - B n (z), E%\z) = A n (z)e 3nz - C n (z). (1.5) 

This is due to the fact that the type II Hermite-Pade polynomials are characterized by a 
Riemann-Hilbert problem, which is directly related to the Riemann-Hilbert problem for 
type I Hermite-Pade polynomials. This relation was first observed in j2Ij and also used in 
[21 E] • See also Section 2 below. So we follow the asymptotic analysis of [U] , and the reader 
is advised to consult a copy of that paper too when reading the proofs in this paper. 

To illustrate the connection between type I and type II we have depicted in Figures ^ 
and [21 the zeros of the type I and type II Hermite-Pade approximants for the case n\ = 
n 2 = n 3 = 60. As it may be apparent from Figures 1 and 2, scaled zeros asymptotically 



TYPE II HERMITE-PADE APPROXIMATIN 



3 




Figure 1 . Zeros of type I quadratic Hermite-Pade polynomials ^60,60,60 

(the 

diamonds on the left) (?60,60,60 (the circles in the middle), and r 60 ,6o,60 (the 
boxes on the right) 



accumulate on specific curves in the complex plane, and the same system of curves is 
relevant for both type I and type II approximants. More precisely, the subsets of zeros on 
the left and on the right in Figures 1 and 2, once scaled, tend to the same limit curves as 
the degree tends to infinity A similar assertion holds true for some of the zeros of g60,60,60 in 
the middle of Figure 1 (those outside of the imaginary axis) and the corresponding subsets 
of zeros of &60,60,60 and C6o,60,60 respectively lying in the left and right half-planes. 



2. RlEMANN-HlLBERT PROBLEMS FOR HERMITE-PADE APPROXIMATION 

In jH] we proved that the type I quadratic Hermite-Pade approximants for the exponential 
function may be characterized as the unique solution of the following Riemann-Hilbert 
problem for a 3 x 3 matrix valued function Y : C \ T, where T is a closed contour in the 
complex plane encircling the origin once in the positive direction, cf. [HI Theorem 5.1], 

(1) Y is analytic in C \ T 

(2) Y satisfies the jump condition 

/ 1 Z -3n-2 e -3nz q\ 

y+(z) = y_(z) o i o , z e r, (2.1) 

\0 Z -3n~2 e 3nz j J 

where Y + (z) and Y-(z) denote the limiting values of Y(z') as z' approaches z G T 
from the inside and outside of T, respectively. 



4 



A.B.J. KUIJLAARS, H. STAHL, W. VAN ASSCHE, AND F. WIELONSKY 



150 




Re z 

FIGURE 2. Zeros of type II quadratic Hermite-Pade polynomials a 60 ,60,60 
(the circles on the left and right outer curves) &60,60,60 (boxes on the right 
inner curve), and c 60j 6o,60 (diamonds on the left inner curve) 



(3) For large z, 

f ( z ' n+1 \ 

z -> oo. (2.2) 




The Riemann-Hilbert problem for Y has a unique solution and Y 21 (z) = P n (z), Y 2 2(z) = 
z~ 3n ~~ 2 Q n (z) for z outside T, Y 22 (z) = z~ 3n ~ 2 E n (z) for z inside T, and Y 23 (z) = R n (z). 

A similar characterization holds for the type II Hermite-Pade approximants. Consider 
the following Riemann-Hilbert problem: find a 3 x 3 matrix valued function X : C \ T — > 
C 3x3 such that 

(1) X is analytic in C \ T. 

(2) X satisfies the jump condition 

/ 1 

X + (z)=X_(z) _^ 3n ~ 2 e- 3n2 1 - z -^-2 e 3nz | ? zeT pa) 
\ 1 

(3) For large z, 

X(z) = (i + O ( - ) ) I I , z^oo. (2.4) 
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Comparing the two Riemann-Hilbert problems it is easy to check that the two solutions 
X and Y are related by 

X(z) = (Y-\z)Y . (2.5) 
The matrix X(z) contains the type II polynomials 1)1 ,4j) . We have the following result. 

Theorem 2.1. Let A n , B n , C n , En , and En^ be as in (|1.4|) and (|1.5j) . where we assume 
that A n is normalized to be a monic polynomial of degree 2n + 2. Then the solution of the 
above Riemann-Hilbert problem for X is unique and is given by 

/ z 3n 2 b ritn , n j r i(?>nz) Q>n,n,n+l(3W'Z) Z 3 ™ 2 C nt n,n+li^ nz )\ 



X(z) 



-3n-2 



BJz) 



Ajz) 



-3n-2 



C n (z) 



[2.6) 



\z 3n 2 b n , n+ltn (3nz) a n>n+ i jn (3nz) z 3n 2 c nin+ltn (3nz)/ 
for z outside T, and 

( —z 3n 2 e^ n ^ n+l {3nz) a n ^ ntn+ i(3nz) — z 3n 2 e| ^ j ^)r^+1 (3n2;)^ 



X(z) 



(2.7) 



\ Z 2e n^n+l,n{3 nz ) 0-n,n+l,n{3nz) Z 3 ™ 2 G^ n+1 ^ n (3nz) J 

for z inside T. In the first rows of (12. 5j) and (|2.7j) we use the Hermite-Pade polynomials 
of indices n,n,n + 1 normalized so that b njnjn+ \(3nz) is a monic polynomial, and in the 
third row we use the Hermite-Pade polynomials of indices n,n + l,n normalized so that 
c n ,n+x,n(3nz) is monic. In addition, we use 

M 



e (1) (z) 



'^ni,n 2 ,n 3 \Z), 



e ni,n2,n 3 { Z ) ~ a ni,n 2 ,n 3 {z)e C ni ,n 2 ,n 3 \Z) ■ (2.8) 



Proof. The entries in the first and third column of ()2.7|) are analytic at z = because of 
the behavior (jl.lj) as z — > 0. So X defined by (|2.6|l and (J2.7|) is analytic on C\T. The jump 
condition (J2.3|) is easy to check with the aid of (jl.5j) and (J2.8|) . The asymptotic condition 
(12. 4|) follows from the degree conditions on the polynomials and the fact that b njritn+ i(3nz), 
A n (z), and c n ^ n+ i^ n (3nz) are monic polynomials. So X is a solution the Riemann-Hilbert 
problem. 

The uniqueness of the solution is a consequence of Liouville's theorem for entire functions. 
Alternatively, it also follows from ()2.5|) and the fact that the solution of the Riemann- 
Hilbert problem for Y is unique, see jHJ Theorem 5.1]. □ 

The fact that the Riemann-Hilbert problems for type I and type II approximants cor- 
respond to each other by applying the inverse transpose is another way to express the 
well-known Mahler relations between type I and type II Hermite-Pade approximants, see 
[TO] or ^21 §2.1]. This simple relation of taking inverse transpose allows one to get the 
asymptotic results for type II approximants directly from the asymptotic analysis that was 
made in [5] for the solution of the corresponding type I Riemann-Hilbert problem. In the 
next section, we introduce several definitions and notations, some taken from jHj and some 
new. Then, we state the results. 
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TZ 



Q 



tz, 



Figure 3. The Riemann surface TZ 

3. ASYMPTOTICS OF TYPE II QUADRATIC HERMITE-PADE FOR THE EXPONENTIAL 

FUNCTION 

3.1. The Riemann surface and related objects. We define TZ as the Riemann surface 
for the function 



z(w) 



1 

w 



1 



1 



w — 1 w + 1 



w(w 2 — 1) 



(3.1) 



The choice of this particular Riemann surface is motivated by the integral expressions for 
the Hermite-Pade approximants to exponentials and the study of their main contributions 
by the classical saddle point method, see [T71 §2.2] or 9, §2.1] for details. 

The rational function (|3.1j) has three inverse mappings. These are the three solutions of 
the cubic equation 



zw 3 — w 2 — zw + — 

o 



0. 



(3.2) 



The Riemann surface TZ consists of three sheets TZp, TZq, and TZr as shown in Figure El 
The bijective mapping tj) : TZ — > C is the inverse of ()3.1j) . We denote its restriction to 
the three sheets by if)p, ipQ, and ipR, respectively. So if>p(z), ?Pq(z), and ipn(z) are the 
three solutions of ()3.2|) . and we assume that the choice of the three sheets is such that 
^p(oo) = —1, iPq(oo) = 0, and ^r(oo) = 1. The Riemann surface has four branch points 
zi = z(w\), z 2 = z{wq), Z3 = z(ws), Z4 = z{w±), which are related to the points w\, W2, w^, 
W4 for which z'(w) = 0. Simple calculations give 



-1/4 



-2k-l 



k 



1,2,3,4, 



(3.3) 
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-3 -2-10 1 

Figure 4. tp- image of the Riemann surface 71 



where ug 
are 



e 2m/8 j g p r j m itive 8th root of unity. The corresponding values of Zk = z(wk) 
2 2 = 3- 1 /VJ, z, = 3-^ull z A = Z-y^ (3.4) 



Z\ — 6 0; 2 4; Z2 — O W245 Z3 — 6 0> 2 4, Z4 — 6 

where co>24 = e 27n//24 is the primitive 24th root of unity. 

The sheets Tip and IZq are glued together along a cut Tp connecting two branch points 
Z\ and z 2 in the left half-plane, and the sheets IZq and Tip are glued together along a cut 
Tp connecting the other two branch points z 3 and Z4, in the right half-plane. Moreover, 
these cuts are chosen so as to satisfy 



7^ [ &Q-ipp)+(s)dseR, 

2vrz J Zl 



zer, 



zeY 



R- 



(3.5) 



(3.6) 



In particular, Tp and Tp are trajectories of the quadratic differentials —{jPq — ipp) 2 {z)dz 2 
and —(ipQ — ipp) 2 (z)dz 2 respectively, see (T3H20]- The -0-images of these curves are shown 
in Figure EJ 

Actually, there are four analytic curves such that (|3.5|) . resp. (13 holds. We denote 
them by Tp, T* P , Tp,i, Tp,2, respectively Tp, T* R , Tp,3, Te,4, see Figure El The particular 
geometry of these curves may be derived from the local and global properties of trajectories 
of quadratic differentials, see [HJ Proposition 2.2]. These curves divide the complex plane 
into different domains. We denote the unbounded domains by -D^p, -Dqo,r, Doo,u and 
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-0.8 -0.6 -0.4 -0.2 



l 1 r 

0.2 0.4 0.6 0.8 



Figure 5. Curves for which either ^ J^ii^Q ~ 4 , p)(s)ds or ^ f* {ipQ 
iPr)(. s ) ds is real. 



-Doo,l, as shown in Figure El We denote by Dp the bounded domain bounded by T P and 
Tp, and by D R , its mirror image with respect to the imaginary axis. 

3.2. Definition of measures and weak convergence of zeros. We now define a mea- 
sure on each of the previous curves. The complex line element ds is defined according to 
the orientation of these curves given in Figure 



Definition 3.1. We define a measure /ip on Tp by 

3 

dnp(s) = —{4> Q - tpp) + (s) ds 
and a measure fiR on Tp by 

3 

dn R {s) = — (V>Q ~ V^)+0) ds 
The measure \ia on Tp U Tr is defined by 

fJ-A = Hp + I^r- 

We define a measure /j^ on TJ, by 



s G T P , 



s G T R . 



(3.7) 

(3.8) 
(3.9) 



d^B(s) = —{4>q - ipp){s) ds 
2m 



s g r* F 



pi 



(3.10) 
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a measure /xc on T* R by 

3 

dfMc(s) = ^—{^Q ~ 4>r)(s) ds s e T* R , (3.11) 

The measures /x B (i> on r^i U Ye, 2 and /x#(2) on r^ i3 U 1^4 are defined as 

3 

dn E a)(s) = — (ip Q - 7p P )(s) ds s e r E)1 u r E)2 , (3.12) 

3 

dii E (2)(s) = —(ip Q -'ip R )(s)ds s G r B>3 U r s ,4- (3-13) 

Z7TX 

The facts that /xp and /xp are probability measures and that fj, E m and /x E (2) are positive 
measures of infinite masses are proven in Theorem 2.4]. So /x^ is a positive measure 
with total mass 2. A similar proof shows that /xp and /xc are positive measures of total 
mass 2, see also the first part of the proof of Proposition 4.8]. 

Denote by VA n i vb„, vc n the normalized zero counting measures of the polynomials A n , 
B n , C n . Furthermore, define the zero counting measures for the remainder functions En^ 
and En^ as 

where the normalization by 2n now corresponds to the degree of approximation and the 
3n + 2 interpolatory zeros of E { n ] and Ei 2) at the origin have been excluded. With all of 
the previous measures, we can state our first result. 

Theorem 3.2. We have 

VA n A ^Ma, "b„ ^ ^Hb, v Cn A (3.14) 

where the convergence is in the sense of weak* convergence of measures, i.e., /x n — > /x i/ 
J fd[x n — > J /d/x /or every bounded continuous function f. Furthermore, we have 

u j?m - ► ^ p (2) — >■ ^(2), (3.15) 

in i7ie sense iTiai 

lim / /(s)di/ = / f(s)dji Ei i)(s), i = 1,2, 



/or ever?/ continuous function f such that f(s) = 0(s 2 ) as s — ► oo. 

3.3. Definition of (/-functions. Next, we introduce the complex logarithmic potentials 
of the measures /xp, /xq, /x^, /xp, and /x^. 

Definition 3.3. We set for S E {P,Q,A,B,C}, 

g s {z) = J]og(z-s)dn 8 (s), z e C \ supp (/x 5 ), (3.16) 
which is defined modulo 27rx. 
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We will also need two more functions, namely 

, , 3 
2 



(ip Q -ipp)(s) ds, 



and 



(iI>q - iI>r){s) ds, 



(3.17) 



(3.18) 



where the paths of integration are in C \ (Tp UT^U {0}). These functions were used in 
[H] and we refer to Lemma 2.7] for their main properties. Note that the logarithmic 
potentials in (J3.16)) have jumps on the support of their measures which can be expressed 
in terms of (fp and <p R , namely we have 



9a+{z) - gA~{z) = -2if P+ (z) = 2ip P -(z), z E T P , 



(3.19) 



9a+(z) - 9a-(z) = -2(p R+ (z) = 2ip R ^(z), z E T 



R, 



(3.20) 



9b+{z) - gn-(z) = -2<p P (z), z e Tp, 



(3.21) 



g c+ (z) - g c -{z) = -2<p R (z), z E T* R . (3.22) 

3.4. Strong asymptotics away from the zeros. As in jH], we will use the function 
V3w 4 + 1 which branches at the four points Wk given in ()3.3|) . We choose as cuts for this 
function the two curves %fjp + {Tp) and iPr+{Tr) (see FigureEJ), and assume that it is positive 
for large positive w. So, in particular we have that V3u> 4 + 1 = —1 for w — 0. 

The following theorem gives the strong asymptotics of the polynomials A n , B n , C n and 
the remainder terms En an d E^ away from their zeros. 

Theorem 3.4. With the functions defined above, we have 

z 2 (^ 2 (z) - l) 2 e^) 



A n (z) 



^q(z) + 1 



1 + 



(3.23) 



uniformly for z in compact subsets of C \ (Tp U T R ), 



B n (z) 



-2) 



n z 2 (ip 2 P (z) -\fe n9B ^ 



l + O ( - 

n 



y/ty P {z) + 1 



for z E C \ (Dp U Tp), 
for z E Dp, 



(3.24) 
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uniformly for z in compact subsets of C \ T* p , and 
. n z 2 (tfj 2 R (z)-l) 2 e n ^ 



C n {z) = { 



;-2)' 



-2)' 



z 2 {^q{z) - 1)V^W 



1 + 
1 + 



Aw*ec\(z$uri) 

/or z G -Dp, 



uniformly for z in compact subsets of C \ T* R . Furthermore we have 

' Z 2 (^ 2 Q (z) - l)2 e n(ffAW-3*) / 
,Z 2 (^ 2 p(z) _ l)2 e n 9B ( 2 ) 



(3.25) 



^0 



0(1 



n 



for z G D^p, 



-2) r 



:-2)' 



i + e>( - 



/orzeC \ (Deep U Dp 

ur P ur Eil ur E , 2 ), 



,z 3n+2 {^ 2 P {z) - l) 2 e - n ^W 



l + O [ - 

n 



/or z G Dp, 



uniformly for z in compact subsets of C \ (T^i U T E2 ), and 



(3.26) 



f z 2 (^qW - l) 2 e n ^)+3z) 



+ 1 



1 + c 



/or z G Dqo.Ki 



( 2) n ^(^)-l) 2 e WflcW A 



V^(z) + 1 



+ I - 



/or z G C \ (A^i? U Dp 

ur R ur £i3 ur Ei4 ), 



;-2) ; 



1 + 



/or z £ D 



R- 



(3.27) 

uniformly for z in compact subsets of C \ (Tg^ U 1^4). 

3.5. Asymptotics near the curves Tp, Tp», Tp, I~p», and r E (i), z = 1,2,3,4. We 
state in the next theorem uniform asymptotics which are valid in neighborhoods of the 
curves where the zeros of the polynomials A n , B n , C n and the error functions En and 
En accumulate, see Theorem 13.21 

Theorem 3.5. Uniformly for z in compact subsets of the region (Tp \ z 2 })UD 00 pUDp 
(this is the shaded region in Figure^), we have 



A n (z) = e n9Aiz) 



z 2 ^ 2 Q {z)-l) 2 



z 2 {^ 2 P {z) - l) 2 e 2n ^) 



y/ty P {z) + 1 



1 + ( - 

n 



where the +sign holds in D^p and the —sign holds in Dp. 



(3.28) 
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Uniformly for z in compact subsets of the region (T R \ {z 3 , Z4}) U U D* r we have 



A n (z) = e n9A(z) 



z\j? Q {z) - If 



I + O ( - 

n 



T + 1 



- o ( - 

n 



where the +sign holds in D^r and the —sign holds in D* R . 
Uniformly for z in compact subsets of C \ D^ p, we have 

B n (z) = _(_2) n e n{9B{z)±Vp{z)) 



z 2 (^(^) - l) 2 e"^) 

" y/3l/}* p (z) + 1 



l + O 



where the +sign holds in Dp and the —sign holds inC \ (A^p U D* p ) . 
Uniformly for z in compact subsets of C \ D^r, we have 

z 2 {i> 2 Q {z) - l) 2 e -^«W 



C n (z) = _(_2) n e n(sc(z)± ^ {z)) 



l + O 



^ 2 (^(^) - l) 2 e»^) 

+ xA^W + 1 



1 + 



where the +sign holds in D* R and the —sign holds in C \ -Doo.r U -D^. 
Uniformly for z in compact subsets of C \ -Dp, we /iawe 



' ^ 2 (^ 2 (g) - / 







* 2 (^) - l) 2 



+ 1 

Uniformly for z in compact subsets of C \ D* R , we have 



\ + o[- 

n 



y/ty%(z) + 1 



O 



(z 2 3^ 2 (z) - I) 2 



3^(2) + 1 



1 + Of - 

n 



(3.29) 



(3.30) 



(3.31) 



(3.32) 



(3.33) 
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Figure 6. The shaded region is where Retpp is negative. It consists of the 
two parts Doo p and Dp, where Dp is bounded by T P and T* p . 

The formulas of Theorem 13.51 are valid on the curves where the zeros accumulate. For 
example, (|3.28j) holds on Tp away from the branch points Z\ and z 2 and (|3.29|) holds on 
Tp away from the branch points z% and z±, that is on the two curves where the zeros of A n 
accumulate. 

It has been proved in Lemma 2.7] that the real part of (pp is negative in D^p U Dp 
and that it is positive in the remaining part of the complex plane, see Figure |H1 For z away 
from Tp so that Keifp(z) < 0, the asymptotic formula ()3.28|) reduces to ()3.23|) . On Tp we 
have Re<£>p(z) = 0, and then the two terms in ()3.28|) are of comparable magnitude. For 
z G Tp, (J3.28j) can be re- written as 

y/Wh-(z) + l V \n)) ^ Q+ (z) + 1 V W;' 

Similar remarks hold for the other formulas. 

3.6. Asymptotics near the branch points. As in jH] the asymptotic formulas near the 
branch points involve the classical Airy function Ai. which is the solution of the differential 
equation y"{z) = zy(z) that satisfies 
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as z — ► oo with |argz| < ir. Because of symmetry, we may restrict our study to the 

asymptotic behavior of A n , B n and En^ near z\. Similar results hold true near the other 
branch points. 

The asymptotics to be given involve a function f\ that is defined in [HI Section 2.7] as 

2/3 



h{z) 



(3.34) 



where we take the branch of <pp(z) which behaves like 

^ p ( z )=c{z-z l fl 2 + 0{{z-z l f/ 2 ) 

as z — y zi, with 0. The 2/3rd power in ()3.34|) is taken so that fi(z) is real and negative 
for z G Tp. Then explicit calculations show that 

f x (z) = c x {z -z x ) + {(z - Zl ) 2 ) (3.35) 

as z — > z\ with 

Cl = f[( Zl ) = 2 1 / 3 3 5/12 e-^ ni . (3.36) 
Theorem 3.6. There is a 5 > 0, such that we have, uniformly for \z — z±\ < 5, 

A n (z) = - Z2 - 1 e - 3 ^e (n+1)(9AW+ ^ W) 



x 



n^hiz) Ai ((n + l) 2/3 /i(^)) (l + O (~) ) 

+n- 1 /%(z) Ai' ((n + l) 2/3 /i(z)) ( 1 + O ( - 



, (3.37) 



-1 p - iw 76 „ /^ e (n+i)(g A (z)+<pp(z)-3z) 



B n ( Z ) = Z- l e- i7r/6 y/7T, 



n^hxiz) Ai (e- 2 «/ 3 (n + l) 2/3 /i(^)) f 1 + f ± 



+n- 1 / 6 / i2 (^)e- 27ri / 3 Ai' (e~ 27ri / 3 (n + l) 2 / 3 / x (z)) (l + O 



(3.38) 



and 



E^\z) = z -l e -™/^ e (n+l)( 9 A(z)+<p P (*)-Zz) 

i^%(z)e- 2 ^ 3 Ai (e^\n + lfl z h{z)) (l + O ^ 

+n- 1 /%( 2 ) Ai' (e 2m / 3 (n + l) 2/3 A(z)) (l + o(- 



(3.39) 



where hi and hi are two analytic functions without zeros in \z—Z\\ < 5, which have explicit 
expressions 

Hz) = [N 2l {z)+ize 3z N 22 {z)) A(z) 1 / 4 , 
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with the branch of the fourth root in fi(z) 1 ^ 4 taken with a cut along Tp, and 

h 2 (z) = (N 21 (z) -ize Zz N 22 {z)} h{zY l l\ 

Here 



2V3#(j) + 1 ' ' " " ^3^* Q (z) + 1 



In Theorem 2.11 of [Oj we used functions hi, h 2 , N 2 \ and N 22 to state asymptotics near 
the branch point z\ for type I Hermite-Pade approximants. The functions hi, h 2 , N 2 i and 
N 22 are the corresponding functions for type II Hermite-Pade approximants. Note that, 
by (J3.19J) . the function g& + fp is analytic near z = z±. 

From the asymptotics near the branch points, one can deduce the behavior of the extreme 
zeros of A n , B n , C n , E n 1} and E n 2) near the branch points. We only state the result for the 
zeros of A n , B n and near z\. Recall that the Airy function Ai has only negative real 
zeros, which we denote by > —i\ > —i 2 > ■ ■ ■ > —i v > ■ ■ ■ . 

Corollary 3.7. Let z^ n , v = 1, . . . , 2n + 2, be the zeros of A n , ordered by increasing 
distance to Z\ . Then for every v G N, we have 

= zi + 2- 1 / 3 3- 5 / 12 e-I-V^- 2 / 3 + O (^j , (3.40) 

as n — > oo. Let z^ n , v = 1, . . . , 2n, be the zeros of B n , ordered by increasing distance to 
Z\ . Then for every i/GN, 



= zi + 2- 1 ' 3 3- s ' 12 e-i"h„n- 2 ' 3 + O (J-Y (3.41) 

as n — > oo. Let z^ , v > 1, be the zeros of En , ordered by increasing distance to z%. 
Then for every v G N, 

z eW = zi - e- 2m/3 - n~ 2/3 + o(-} 

f>(zi) n +U {n) 

= zi + 2- 1 / 3 3- 5 / 12 eI™^- 2 / 3 + O , (3.42) 



as n — ► oo. 



We note that the previous formulas are exactly the same as those obtained for the zeros 
of the Hermite-Pade approximants of type I, P n , Q n , and E n near zi in Corollary 2.12]. 
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4. Proofs of the asymptotic formulas 

4.1. Transformations of the Riemann-Hilbert problem. In [9J, the analysis of the 
Riemann-Hilbert problem for the matrix Y was made through a sequence of explicit trans- 
formations 

Y i-> U ^ S, (4.1) 
leading to a matrix valued function S that satisfies 

S(z) = I + o(i) (4.2) 

uniformly for z G C\ Eg, where is a system of contours shown in Figure 13], see [HI 
Theorem 6.4]. 

The starting point of the asymptotic analysis is the Riemann-Hilbert problem for Y 
on the contour T as given in the beginning of Section 2. We choose T so that it passes 
through Tp and Tr, and the remaining parts are in D^u an d -Doo.l, see Figure El The 
transformation Y i— > U is given by 

= L- n_1 y(z)diag ( e " (n+1)9p(z) , e (n+1)sA(z) , e -( n+1 )^)) L n+1 (4.3) 

where L = diag (e £//3 , e~ 2i ^ 3 , e^ 3 ) with £ = log 2 — ni, see Section 5.2]. This transforma- 
tion has the effect of normalizing the Riemann-Hilbert problem at infinity (i.e., U(z) I 
as z — > oo) and to introduce oscillatory jumps on the contours Tp and Tp. 

The transformation U *—>■ T consists of opening of lenses around Tp and Tp. We will 
not give all the details about this transformation here, but we refer to Section 5.3]. 
The result of the transformation is that T is the solution of a Riemann-Hilbert problem 
on a complicated system of contours, but the jump matrices on all contours tend to the 
identity matrix as n —>■ oo, with the exception of the jump matrices on Tp and Tp, which 
are independent of n. 

For the final transformation T ^ S we surround each of the branch points Zj by a small 
disk, and we define S = T(M^)~ 1 inside the disk around Zj, and S = TN' 1 outside 
of these four disks. The matrix- valued functions N and j = 1,2,3,4 are explicitly 

constructed parametrices that resemble T as n gets large. The result is that S = l+0(l/n) 
uniformly 

Since X = Y~ f , we get from (|4.1|) the corresponding transformations for X 

X^U-*^ T-* i-> S'\ (4.4) 

and from ()4.2j) we also get S~ l = I + 0(1/ n) uniformly on C \ S5. The proofs of the 
theorems then follow in the same way as the proofs in |9, Section 7]. The role of the 
exponential prefactors is now played by N := iV~*, and so to get the explicit formulas of 
Theorem 3.4 we need to evaluate the entries of N. 

4.2. Calculation of N. We have that N := N~ l solves the following Riemann-Hilbert 
problem (compare [HI §6.1]), 

(1) iV is defined and analytic in C \ (Tp U Tp). 
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(2) N has jumps on T P and T R given by 

N+(z) = N-(z) 

and 

N + (z)=N_(z) 

z -iJ* 



z^e- 32 








— ze~ 





zer. 



z g r 



(4.5) 



(4.6) 



(3) jV(z) =/ + C(i) asz-^ oo. 
Proposition 4.1. v4 solution of the Riemann-Hilbert problem for N is given by 

(HM*)) li(M*)y 

N(z) = F 2 (Vp(*)) F 2 (if> Q (z)) Fj{^ r {z)) 
\F^ P { Z )) F 3 (i> Q (z)) F^ R {z)) t 



where 
F(w 



—w(w — l)G(w) 
V3w 4 + 1 : 



Fo(w) 



{w 2 - l)G(to) 
3V3W 4 + 1 



F*(w) 



w(w + l)G(w) 



(4.7) 



(4.8) 



V3w 4 + 1 

u>i£/i \/3w 4 + 1 defined and analytic in C\(ijjp + (T P )Uip R+ (T R )) , and such that it is positive 
for large positive w. The function G is defined by 



G{w) 



( m + l)(2m-l) 



«i 2 -l/3 



g U^ — 1 



(w-l)(2m + l) 

V. w 



for w G ip(7Zp), 

for w G VC^-q), 
/orw G ip(U R ). 



(4.9) 



Proof. Rather than taking the inverse transpose of the solution found in jHJ Proposition 
6.1], we prefer to repeat the argument given there in the present situation. One checks 
easily that the only difference lies in replacing the jumps in Equation (6.14)] by their 
inverse. Hence the function G used in the proof of 0| Proposition 6.1] should be replaced 
with G = G^ 1 . This leads to the asserted solution for N. □ 

Corollary 4.2. The entries in the second row of N can also be expressed as 



N 21 {z) 



N 22 (z) 



N 23 (z) 



2y/ty*(z) + 1 

z\r Q {z) - l) 2 . 



(0 



R\ 



3l^(*) + 1 
iy e gR(z) 



2y/3^ R (z) + 1 



(4.10) 



(4.11) 



(4.12) 
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Proof. From ()4.8|) and the definitions of the functions F 2 and G in Proposition 6.1], we 
derive that 

~ ( W 2 ~ 1) (w 2 - l) 2 

p (iu) = — = . 

V ; 3V3«; 4 + lG(w) (3w 4 + l)F 2 (w) 
Now, replacing w with ipp(z), iPq(z), ipp(z), respectively, and using the expressions of 
N 21 (z) = F 2 {il> P (z)), N 22 (z) = F 2 {il> Q (z)), N 23 (z) = F 2 (^ R (z)) that are given in Equa- 
tion (6.20)], and also noting that gA = 9p + 9r, we arrive at (|4.1U|) - ()4.12|) . □ 

4.3. Proof of Theorem EIH We get from (JOJ) and (JO}, and the fact that U = T = SN 
away from Tp and Tp, that 

X =y~* = L- n - 1 (/ + C(l/n))iVdiag (e~ (n+1)9P , e {n+l)9A , e - {n+l)9R ) L n+l (4.13) 

uniformly in compact subsets of C \ (Tp U T R ). 
From (ETT3J) and (|2~6>(|2~7)l we get that 



X 22 (z) = N 22 (z)e^ n+1 ~> 9 ^ (l + (^j , 



uniformly for z in compact subsets of C \ (Tp U Tp), which by (|4.11|) leads to (J3.23j) . 
In the same way, we get from 1)4.13)1 and (J2.6)) that 

B n (z) = X 21 (z)z 3n+2 = N 21 (z)(-2) n+1 z 3n+2 e^ n+l)9p{z) (l + O (-Vj , 

uniformly for z in compact subsets of the exterior of T. Using ()4.10j) and the fact that 
9b(z) + gp(z) = 3 log z for z G C \ Dp, we find the formula f)3.24|) for z in the exterior of 
T. For z inside T, we have by ()4.13|) and ()2.7j) 

B n (z) = A n (z)e' 3nz - E£\z) = X 22 (z)e~ 3nz + X 2l (z)z 3n+2 

= N 22 (z)e^^e-^ (l + O ) {m 

+ N 21 (z)(-2) n+1 z 3n+2 e- {n+l)9p{z) (l + f ij V 

The first term in ()4.14|) dominates for z G Dp, and this leads to the second asymptotic 
formula in ()3.24|) . since gA(z)—3z = gp,+£ for z G Dp. The second term in ()4.14|) dominates 
in the remaining part inside T and this again leads to the first asymptotic formula in ()3.24|) . 
A more detailed analysis for z near Tp and Tp reveals that ()3.24|) is uniformly valid near 
Tp and Tp as well. 

The proofs of lET25|) . and (ET27I) are similar. □ 

4.4. Proof of Theorem 13.21 The limits for the counting measures z/^ n , VB n , and vc n 
follow directly from the strong convergence results in Theorem 13 . 41 and the unicity theorem 
for logarithmic potentials, see e.g. fSJ Theorem II.2.1]. 

To establish the limits for v „(i) and v„{i) we follow the proof of 19, Theorem 2.5] given 
in Section 7.2 of [9 , which established the limit for vp n as n — > 00. In following this 
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proof we have to make some obvious modifications to see that the proof applies to the 
present situation as well. The only modification that may not be immediately obvious is 
that in (HI Equation (7.13)] we used an integral representation for E n to conclude that 
- log(z~ 3n ~ 2 E n (z))\ z=0 is bounded from below as n — ► oo. We do not have an integral 

representation for En and E n 2 \ but the fact that - \og( z~ 3n ~ 2 En (z)) is bounded 

n 2 = 

from below as n — > oo for j = 1,2, follows easily from the asymptotic formulas 1)3.26)1 and 

flUZI). □ 

4.5. Proof of Theorem 13. 5L This is similar to the proofs of Theorems 2.9 and 2.10 in 
|Hl Sections 7.3 and 7.4]. We omit the details. □ 

4.6. Proof of Theorem 13.61 The proof is similar to the proof of Theorem 2.11 in [§J 
Section 7.5]. Applying the inverse transform to [9, Equation (7.19)] leads to an expression 
of X(z) involving the matrices N and (S~ l Y = I + 0(l/n). Multiplying that expression by 
the unit vector (l 0)* as in the proof of Theorem 2.11] leads after some calculations 

to the asymptotic formula ()3.39)) for E n , while multiplying by the unit vector (0 1 0)* 
leads to the formula (j3.37j) for A n . In order to get the asymptotics ()3.38|) for B n , one has 
to choose z outside of the contour V for the Riemann-Hilbert problem for X. Then one 
performs similar calculations in the regions / 1 _1 (J) or f~ l (II), as they were defined while 
constructing the parametrices near the branch points in §6.2], see also [3 Figure 12]. 
We note that the following relation between the Airy function and its derivative is used in 
the above computations, 

gWr/6 

Ai(s) Ai'(ws) - u 2 Ai(us) Ai'(s) = — - — , 
where u = e 2ni/3 , see e.g. jEH Ex. 8.1, p. 416]. □ 

4.7. Proof of Corollary 13.71 The formulas for the zeros of A n , B n , and En near z\ 
follow from the asymptotics near the branch points given in Theorem 13.61 We refer to the 
proof of Corollary 2.12] given in [HI Section 7.6] for details. □ 
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